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ABSTRACT 


The persistence of the free expansion of the detonation products of a 


high altitude detonation is determined by the sum of the stresses of the ,en- 


vironment. Below 200 k m the classical dynamic friction with the air couples 


the detonation products to a very much larger m a s s of air. Above this alti- 


. 
. 


tude more uncertain plasma phenomena. of magnetic and electrostatic shocks 


must be considered to determine this coupling. 
If a cubic scale height of air 


1 


i s raised to a temperature such that its sound speed is greater than escape 


velocity, it jets upward, and the ionized fraction is stopped and mixed un- 


stably with the earth's magnetic field. The subsequent expansion along the 


field lines deposits the major fraction of the debris a t the opposite conjugate 


point. he P decay in transit on the various L surfaces reached by. unstable 


mixing determines the high energy electron injection into trapped orbits. For 


the Starfish event (1.4 megatons at 400 km) theory predicts that the detonation 


products and a cubic scale height of air expands upwards until stopped by the 


earth's field at, 1000 km. In addition one, expects the deposition of the debris 


in the atmosphere at the southern conjugate point at 250 km.altitude dis- 


persed over 300 k m due to magnetic field mixing. Approximately 5% P injec- 


. 
. 
. 
. 


tion into trapped orbits of L > 1.3 is expected. 
. 
. 
. . 
- iii - 
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The Phenomenology of the Mass Motion .of 


a High.Altitude Nuclear ~ x ~ l o s i o n " 


Stirling. A. Colgate 


Lawrence Radiation Laboratory, University of California 


Livermore, California 


May 13, 1963 


The initial mass-motion of a nuclear explosion is, of course, governed 


by the hydrodynamics .of th.e free expansion of the detonation products; how- 


ever, the subsequent persistence of this motion is determined by the sum of 


the stresses of the environment. At a high enough altitude the expanding 


detonation products may encounter less than.an equivalent mass of a i r and. a 


free expansion will take place until s'ome additional stress,. a s , f.or instance, 


magnetic field, becomes dominant. 


ATMOSPHERE DOMINATED EXPANSION 


The pertinent question that separates the behavior of the free expansion 


from that governed Ly environmental density is: at what altitude will the 


integral m a s s of air above the explosion that ' collides ' with the detonation 


products be greater than the m a s s of the products ? 
Since the atmospheric 


density near the earth varies exponentially with height, the supersonic 


momentum- conserving expansion for a snowplow-type shock has a velocity 


distribution 


- 


:k Work done under the auspices of the U. S. Atomic Energy Commission. 
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where 


R = radius of expansion 


M = m a s s of detonation products 


uO = initial mean expansion velocity 


8 = angle with respect to the vertical 


and 


~ ( 0 ) 
= p0e -(r cos 8)/h 


where 


h = scale height of the atmosphere 


p o = atmospheric density (in g/cc) at the detonation point. 


The shape of this shock front at some later ti1.1.1e t is 


If we look ~ L L Lllis sliapc aftcr a time long enough s n t h a t R(f3 = 0) >> h 


and concern ourselves with cases where the total grams per cmL of air is 


large cnmpared t o the detonation product m a s s at a radius large compared to 


h, namely, 


then the radial distribution becomes, Fig. 1 


.M cos (3 
R ( 0 , ,t ) - u 1. ------..- 
.. 
h 
0 
3 
for 1 - COG e /; - 


8.rrh P o 
R(O) ' 


This distribution is based upon the assumption that the detonatiorl 


products ' sweep up1 the air ahead in a lsnowplowl fashion, i. e., a strong 
- 


shock with the adiabatic Y -, 1. If, on the other hand, the detonation product's 


could interpenetrate the ambient air with negligible interactibn, then very 


obviously the free expansion approximation would apply. 
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DETONATION PRODUCT-AIR INTERACTION 


Three interaction mechanisms between the detonation products and the 


air will be c-onsidered: ( 1). the Coulomb c.ollision dynamic friction (both ion- 


ion and ion- electron), (2) magnetic shocks, and ( 3 ) electrostatic shocks. 


Each of these interactions depends upon the initial fractional ionization 


of the air. Two processes for ionizing.the air a r e evident: first, the pulse of 


X rays will give an initial partial ionizati'on, and second, each of the above 


mechanisms, of interaction will lead to electron heating. and subsequent 


ionization. The energy required for complete ionization represents a small 


fractional coupling for any of the above three mechanisms since the energy to 


be exchanged is very large compared to ionization potentials. In those cases 


where the initial condition for strong interaction requires a modest fractional 


ionization the X-ray pulse will satisfy this condition. 


F r o m Eq. (5), the effective path in which interaction must take place is 


approximately h. 
The flux of 1-kev X rays from a nominal 1 megaton of 


1-kev X rays a t a distance h, assuming small attenuation, is 


where 


WB = 1 megaton 


W = 1 kev 
X 


h = scale height = 50 km. 


For a mean absorption coefficient wf 1 eke v X rays in nitrogcn and 


d 
-20 
2 
oxygen of u = 5 X 10 
cm the available energy per atom becomes 


Wa = +uW 
= 4 ev. 
X 
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The 'energy invested per ion pair is. usually 32 ev, and considering 


that the recombination time at the altitudes in question is longer than the 


delay for the arrival of the detonation products At x h/u 
= 0.015 sec, then 
0 


a fractional ionization of 12% or greater should exist a s a minimum initiation 


level for any of the str.ong interaction mechanisms. 


COLLISIONAL DYNAMIC FRICTION 


The'dynamical friction of a test ion s'lowing down in ;i, plasma has h 
~ 
~ 
n 


. 
. 
treated in detail by Kranzer [~ranze?, 1961 1, who c o n f i ~ ~ m s 
'tlie slvwing dowll 


approximation of Spitzer [ ~ ~ i t z e r , 
1962 1. 
The dynamic .friction on a test ion 


can be separated into two parts - that due to collisions with ions and that due 


to c.ollisions with electrons. The electr.on temperature is a contr.olling fac- 


tor for the electron dynamic friction, because it governs both the thermal- ' 


ization time and the charge state of the ions. 
(The usual pickup and loss 


dynamic. friction for heavy ions passing. through neutral matter is smaller 


than the free electron friction for electron temperatures > 5 cv.) 


The critical region for determining interaction is that region a scale 


height from the detonation where the products arrive essentially noninteracted 


with a velocity u and spread in time h/u 
secondo. The .nlaxirnurn heating 
0 
0 


rate then becomes 


m e 
2 
2NZ 
M-u h 


dW ' 
. 
L 0 


- = 
z 
dt 
erg/cc. sec. 
' 
(7) , 
m e 1h(kTe)3h/~e4neZ2 gnA 
- 


where NZ in the m~rnher density of thcdetuliation products, M 
its atomic 
i 
: 
. 
z 


C 
m a s s , Z its charge, and n the plasma electron density. 
e 


The loss rate of energy from the plasma wi1.l be entirely by line emis- 


sion (brems s trahlung being negligible). An .approximation to this 10s s rate 
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can be made according to Post [post, 
- 
19611 and an independent evaluation of 


the Model C Stellarator impurity radiation loss rate st off man et al., 19621 


can be made by assuming that each atom which is not stripped below three 


electrons has an allowed level 10 ev above the ground state with an excitation 
-16 
.2 
cross section u = 10 
cm . This gives rise to a loss rate 


- - no z ( k ~ / 2 m ) l / ~ 
lo-15 ev/cc- sec. 


An.approximation to the ionization state Z can again be.made from the 


work of Post a s 


and 


Letting.AZ = 50 and 


then, equating heating and loss rates, one obtains 


or 
11 
10 
11 
T = 6 X 10 
/no ev, for 1.2 X LO 
< n < l.2.X LO .. 
- 0 . - 
. 
(10) 


For densities gredter than l.2.X 10''/cc the approximations (8) for the 


resonant radiation loss rate and (9) for the effective Z of the bomb material 


break down. Since the resonant-state energy is 10 ev, temperatures below 


5 ev give ris e to an .exponentially reduced radiated power, effectively es - 


tablishing .a lower .limit to the electron .temperatur.e of -5 ev. In addition, 


the pickup and loss cross sections of the detonation products traversing 
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8 
neutral matter at 3 X ,10 cm/s ec estab1ish.a.minimum effective Z of 


approximately 2. 


4 


At the critical air density, below which there is no interaction, the 


dynamic friction must b e sufficiently large to slow down the detonation 


products ions in a time t = h/uo. 
Therefore 


1 12 


- M~ me 
(kTe)3 I2 
h/uo - ;;;- 


, 
, - 
4 
2 
s ec. 


e 
r e n .Z -1nA 


0 


2 
8 
UsingEq. (9) for n e and Z , and uo = 3 X 10 cm/sec and h = 5 
6 


X 1 
c m , then 


M2 me 


/2 (k)3/2 


1.6 X lom2 = ;;;_ 
4 
where dn A = 15. 
e r e n I n A 
0 


Then 


n = 1% ~ o ' l / c c . 
0 


I. i 
l h e nuclear Coulomb scattering cxiterin,n f ~ ; r 
detonation producto air 


interaction requires .a still higher density: 


6 x 
Z12.zf 


n h m = 
2 


U 
c 


air 
8 
-19 
2 
For u 
= 3 . X 10 cm/sec, 
cr 
= 6 x 10 
c m 
0 
C 


Thus the electron dynamic fiiction is copsiderably larger than the nuclear 


coulomb scattering. 
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MAGNETIC SHOCK IN.TERACTION 


The strong interaction by means of a magnetic shock wave depends 


upon the momentum stress .of the detonation products compres.sing.and moving 


the ambient ionized air. If we first neglect the presence of .the ionized air 


and ask how the detonation pi-oducts interact with the field alone, then the 


simple c.oncept of the ionized detonation products blowing .a diamagnetic hole 


in the .magnetic field is valid. In general, the hole will expand until the work 
. 
. 


done against the' field equals the kinetic energy of the detonation products. 


That is, 


3 
2 
22 
4 / 3 ~ R : B /BT 
= 4 X 10 
erg 
8 
R = 10 cm.for B = 0.3 gauss. 
(13) 


The thickness of the boundary X between the diamagnetic sphere of expanding 


ionized detonation products and the magnetic field 'is described by the ' M 


theory ' of Ros enbluth '[Ros enbluth, 


2 
2 1/2 
X = (mc /8rr ne ) 


where m = electron mass, and 


n = number density of electrons associated with detonation 


products at the boundary. 


Since the field can .move at the velocity -of light, and the detonation 


products move very much more s.lowly, there will be a low density front 


where the momentum. s t r e s s equals the field pressure. This is the condition, 


for which the above thickness i's derived. The.physica1 description of such a 


layer is that the ions penetrate a distance X into the field with negligible 


interaction, but the electrons associated with the ions cannot penetrate the 


field a s far because of their smaller Larrnor radius. The resulting charge- 


separation electric field decelerates the ions and causes.a transverse drift 
.." - 
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of the electrons. This electron drift gives rise to a current which is.just 


the diamagnetic current bounding the field from the diamagnetic hole. 


For the case of expanding detonation products a t a distance h spread 


>. 


out in time h/u 
the maximum magnetic field interaction distance X becomes; 
0 ' 


2 
X = 2500 c m , where n M u 
/2 = B ~ / ~ T 
z 0 


4 
n = 10 for B = 0.5 gauss 


This then becomes the maximum length that the detonation products can 


penetrate the plasma unless an instability. takes place which causes a break- 


down of the layer. There a r e two types of instability that can allow the 


plasma to penetrate the field: one is the hydromagnetic Taylor instability, 


which depends upon the deceleration of heavy fluid by a light fluid, and the 


second is an electron velocity phase-space instability in which the electron 


drift stream is unstable to the generation of plasma oscillations. 


The hydromagnetic instability obviously takes place whenever the 


detonation products density is g;r eater than .the field-plus -r eoidual-plasina 


density. As has been shown experimentally, [ ~ i c k i n s o n et al., 1962 1, one 


would expect this instability for the conditions of free expansion whenever 


the field stress decelerates the detonation products. For the case of the 


earth Is field this occurs at R = 1000 krll and so does not affect the initial 


question of coupling to the air. (It will of course, affect the final plasma be- . 


havior that we will discuss later.) The case of coupling to a 1arger.mass of 
. . 
air should therefore be Taylor stable.. 


The sheath inctability [Dunelrla~l, Lu be published; Hartman and Sloan, 


to be published; Shkarofsky, to be published] takes place when the electron 


Larmor radius becomes equal to the Debye length by the generation of plasma 


oscillations which allow the plasma to interpenetrate the field. However, this 


-9- 
UCRL-7224 Rev. II 


interpenetration can take place only at a rate given a s . a maximum by Bohm 


diffusion DB = ckT / e ~ . Since this maximum diffusion corresponds to a 
e 


random walk step of a Larmor radius in a Larmor period (ions and electrons), 


then the effective diffusion velocity must always be small compared to the 


directed s t r eam velocity. Consequently, regardless of the instability a strong 


momentum interaction.between the driving detonation products and the driven 
i 


air plasma will have taken place. 
The instability merely affects a major un- 


settled problem concerning the shock structure [ ~ a r d n e r et al., 1958; Adlam 


and Allan, 1958; Colgate, 1959; Auer, Hurwitz, and Kolb, 196 1; Morawetz, 


19621 namely, the immediate division of internal energy among electrons, 


ions, .and ion oscillations. 
For complete coupling, however, the strength of 


this interaction requires that either the s t r e s s is spread out by instability over 


a distance large compared to:X, or that the ambient magnetic field i s com- 


pressed to-a value so that its pressure equals the time rate of change of the 


detonation products momentum. 
That is, 


3 
Wt /8n h 
= B 2 
max /8 7-r 


B 
= 20 gauss. 
max 


This is not too large a compression, but a shock structure of such a high 


~ l f v e / n 
Mach number has not been derived on a self-consistent basis. 


NEUTRAL STRESS MAGNETIC SHOCK 


A more conservative picture of the magnetic shock interaction can be 


formed based upon a time-dependent ionization model similar to.that formu- 


lated by Gerry, Kantrowitz, and Petschek err^, Kantrowitz, and Petschek, 


19631. In this model the time-dependent ionization of the neutral fraction oc- 


curring within the moving combined shocked plasma and magnetic field r e - 
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sults in the primary momentum stress (Fig. 2). The effective shock thick- 


ness becomes the ionization mean free path, which may be longer than the 


instability wavelengths discussed previously. 
As a consequence the shorter 


wavelength instabilities can be neglected, and the ionization mean free path 


in the shocked plasma becomes the critical parameter that determines 


coupling. The ionization of the neutral air occurs by: 


1. The X-ray pulse (Eq. (6) - the electrons from which may or may 


not have reached equilibrium by the time of arrival of the shock); 


2. 
The electron ionization; 


3. The charged-detonation-products ions and shocked-air ions. 


Even assuming equilibrium following .the X-ray pulse, this .fractional 


ionization is small compared to unity and so the primary coupling must be 


determined by 2 and 3. 
For electron.te'mperatures above 5 ev, th,e ionization 
, 
. . 


... 
, 
' 
-7 
rate of a i r is approximately constant [ 
a>$-Ma-ss-ey, 
19491 ( q v ii 3 


. - 8 ' 
3 
X 10. 
cm /sec), and s o electron ionization should be indepe~ldent of shock 


velocity for air densities where magnetic shocks a r e important. The ioniza- 


tion by partially stripped ions, on the other hand, is sensitively dependent 
, 
, 
, 


. . . 
upon velocity. Choosing the initial detonation product velocity as the cr.i.tica1 


velocity for coupling then the equilibrium char,ge state must be determined 


I 
. 
, . 
. 
for the ionization rate'. Knipp and Teller [ ~ n i p p 
and Teller, 194 1 ] haye 


I 
. 
: 
" 
. . 
estimated the effective charge of a high-Z atom slowing down in a ne,utral gas 
. . 


and have shown that capture and 10s s' a r e in eiquiljbrium when u = v /1.5, 
0 
e 
8 
where v is the bound orbital electron velocity. 
Choosing u = 3 X 10 cm/sec, 
- 


e 
0 
8 
then the ionization energy corresponding to v = 4.5 X 10 cm/sec is 80 ev or 
e 


= 3.5. 
'eff 


The cross section for nitrogen ionization by an ion of charge Zeff and 
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The detonation product density at a distance 'h is 
. . 


. 
. 
3 
6 


N~ = M X 6 X 1oZ3/x 4/3nh. 
For M = 10 gram and A = 20, 


Therefore the mean free path A of a neutral air atom in the flux of detonation 


products becomes 


6 
. A = NU = 5 X 10 cm. 


The electron density such that electron ionization will have contributed 


equally to the ionization of the air within the time of the passage of the detona- 


tion products is: 


In summary, a strong interaction between detonation products and' the 


combined air and magnetic field is predicted for an ionization-type hydro- 


magnetic shock where the shock thickness due to either ions or electron 


ionization is large' and 'comparable to the scale height of the atmosphere. 


During,the expansion phase such a large shock thickness should be stable un- 


til the magnetic field stress reverses the expansion, 


Since the ionization magnetic shock by debris ions becomes the dominant 


shock mechanism above'.350 km, it is pertinent to demonstrate that the high 


debris charge state of 3.5 may have had a chance to reach.and maintain 


, 
, 
equilibrium before having expanded to a s cale height in radius. It i s evident 


that a free expansion in a vacuum' will permit the detonation. products to cool 


by adiabatic expansion and radiation leading. to recombination and a lowering 


of charge state. However, because of collisions with air molecules the 


detonation products will b e heated,, but this heating is small ' compar ed to 


' blackbody radiation (even at 1 ev). ' Consequently, it is only after the detonation 
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products-have expanded to the 
where they a r e optically thin, yet aq air- 


atom collision.range thick, that the temperature will r i s e and the charge 


state ,become high. 
Naturally, if the detonation .products a r e a collision 


range thick, there is more than ample time for charge equilibrium to be 


reached. 


The range from.Eq. (1 1) becomes 


where 


ND = detonation product atom density 


R = radius of -detona.tion products. 


Solving ,for R analogously to Eq. ( 16) gives 
12 
R = 1, km, ND = 5 X 10 
/cc. 
- 


Using Eq. (8) for the resonant line emission power loss (an upper limit 


because of the available states and self -absorption) ,and assuming 


n 
= 3.5 n gives a power loss rate 
C 
0 
2 
= 10" ergs/cm s e c . 
',ad 


The a.ir density p 
such that the collisional heating e'quals the power radiated 
a 


At lower air densities the balan'ce will occur at a larger detonation product 


radius since the heating rate .decreases a s 1 / ~ 
for R > air atom range; 
D 


L 
whereas. the radiation decreases a s 1 / ~ . Consequently, over the air 
D : 


-11 
density range of interest ( L O 
2 p a 2 1 0 - l ~ ) 
the detonation-product charge 
. 
, 
. 
I 


state should be maintained at an equilibrium value of Z z 3.5 until slowed 


down, 
. . 
I .. 
I . . . . 
... : 
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ELECTROSTATIC SHOCK INTERACT ION 


i 


' 
. The electrostatic shock depends upon the growth ofilarge-amplitude 


plasma os ciilations, and/or ion sound waves due 'to the two- stream interaction. 


Buneman [ ~ u n e m a n , , 19591 has shown that the electron-ion two-stream in- 


stability grows to saturation in roughly 30 plasma periods. Saturation occurs 


when the electron stream kinetic energy is randomized by nonlinear effects 


to a quas i-thermal distribution. When the r esulting electron thermgl velocity 


(Veth ) equals the relative stream velocity,. linear-unstable growth stops. The 


ion-wave instability a s shown in detail by Stringer [stringer, 196 1 ] grows .pro- 


vided the ion temperature is small compared to the electron temperature 


, . 
and provided the electron thermal energy is l'ess than the ion' stream energy. 


On the other hand, the maximum ion wave amplitude (measured in energy) is 


limited in linear theory to the electron temperature. 
The possibility of 


further .nonlinear instability growth has been investigated numerically by 


,Hartman and Colgate a art man and .Colgate., 1962 1, and although large 


amplitude growth has been confirmed for ele,ctron .ion mass ratios up.to ,100, 


so+e doubt exists for m a s s ratios a s great a s lo5 for the present case. 


_ 
a 
' 
/ 
Accepting the more conservative result of the saturation of the electron 


instability a t v 
= u 
this gives a two-stream dynamic friction.that will 
eth 
0' 
2 
maintain the .electron temperatures at' m u 
/2 ,-. 50 ev. 
This ma'gnitude of 
0 


the friction is then,limited by the 'plasma radiation loss rate a t 50 ev. 


Equating ,this loss rate to the detonation products energy input rate determines 


a minimum density for .interaction. From Eqs. (8) and (9), 


UCRL-7224 Rev. I1 
- 14- 


where 


1 /2 
(hv) = 10 ev, ( k ~ / 2 m ) 
= uo 


8 
= 3 X 10 cm/sec, T = 50 ev, and 


Then 
10 
0 = 1.4X 10 
/cc. 


This is a stronger interaction than the electron Coulomb dynamic friction, 


but certainly weaken than the magnetic shock. 


ENERGY CONSERVING INTERACTION AND EXPANSION 


The magnetic shock ensures strong interaction essentially independent 


of density, but the limit of free expansion behavior is determined by the radial 
6 
distribution of Eq. (5);: hooki in^ M = 10 grams, and h = 50 km, then pmin 


becomes s 


Above this altitude free expansion takes place. 
Below this altitude the com- 


bination of detonation products and heated air will expand upward at a velocity 


given by Eq. (1) provided the resonant-line' radiation cooling of the a i r behind 


the shock justifies the approximation of a snowplow shock, that i s , that the 


internal energy is radiated away a s fast a s it forms. However, for densities 
7 
'. 
l e s s than,the limitsof Eq. (lo), that is, lo1' < n c 10 /cc, the resonant 
0 
radiation 10s s very rapidly becomes negligible, because air atoms a r e stripped 


of electrons below the lithium 3-electron configuration, and only levels of 


I 


large energy gap and small cross section remain. The radiation power loss 
. 


rate f r o m a i r decreases by two orders of magnitude as the temperature in- 


creases from 20 to 100 ev [ ~ o s t , 19611, 
consequently, in this intermediate 
- 


15- 
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density region the subsequent expansion of the high-temperature air plasma 


- 


'controls the.hydrodynamics more than the directed momentum of the detonation 


products expans ion. 


In order for this 'bubble' of high temperature plasma.to expand up- 


wards through the remaining atmosphere, it 'must comprise a volume of 


2 
roughly a cubic scale height, otherwise the mass/cm ,above will be greater 


than .that in the bubblef, thereby limiting upward motion. 
The expansion 


velocity upward is then sound speed within the bubble. Assuming a specific 


heat per electron .of 3 to include the heat of ionization and the radiation during 


ionization and a volume from Eq.. (5), the expansion velocity u 
becomes 
ex 
3 
24np h (ueXl2/2 = Wt 


u . - 
ex - (W t /12np h3) 
cm/s ec. 


The bubble will cool by both radiation and expansion, but in general 


once the temperature has run away above the resonant line-emission range, 


the adiabatic. cooling will dominate. As pointed out earlier, the subsequent 


expansion will be limited by the stress of the magnetic field at roughly 1000 


km for a megaton. Assuming a V = 5/3 adiabat, the final temperature T be- 
f 
comes 


whcrc 
n 
= 5.2, n for T. > 50 ev. 
e 
0 
1 


Then 


3 
For R = 10 km, h = 50 krn, W = 1 megaton, and 
t 
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When the bubble material. is stopped by the magnetic field, one would expect 


some heating of the electrons from the inelastic energy of the ions. 
The 


'fractional cooling to be expected within the time of interaction between the 


bubble material and magnetic field (assuming T lies within the resonance 
e 
radiation range) by adiabatic cooling i s by Eq. (8): 


This fraction is small (- 
for the largest no (- 10 l o / ~ ~ ) , 
and so the 


radiative cooling of the plasma after expansion is negligible. So also is the 


heating for T - 
. 
> T 
. T max is that electron .temperature such that the 
max 


electrons can make one thermal exchange time with themselves within thc 


expans ion time. 


T = u.3 n /u 
ev. 
0 
ex 


When the adiabatic temperature of Eq. (22) exceeds the heating temperature 
10 
given by Eq. (24), then no further change should be expected, At n 
= 10 
/cc 
0 


the electron temperature should stabilize at roughly 10 ev from both effects. 


MAGNETIC FIELD LIMITED EXPANSION AND FINAL DE.POSIT ION 


When the plasma of a i r and detonation products is stopped by the stress 
\ 


of the magnetic field at a distance R, we expect a Rayleigh-T.aylor unstable 


mixing ,between the two because- the density of the incident p 1 a s m a . i ~ greater 


than.the ambient density (Fig. 3). The linear growth of a perturbation'of 


amplitude A, .is [ ~ i c k i n s o n et al., 1962 1: 


' '/2 t ] 
A = A, exp'[g./2n h ) 
(25) 


where g is the acceleration and A the wavelength. 
For a linear deceleration 


taking place in a distance R, the growth becomes: 


1 /2 
A = A, exp [ ~ / n 
A ] 
. 


Only the nonlinear phase of instability growth (where .A - 
> ~ / 2 ) 
will irreversibly 


mix the- plasma with the fie1.d. Conscqucntly, if we wish to lcnow the i r r e - 


ver sible -mixing of the plasma onto magnetic flux surfaces other than the one 


defined by the explosion, we must ask what is.the largest wavelength that r e - 


quires the fewest generations of growth to reach an amplitude greater than 


~ / 2 . Then, because of the.faster growth rate of smaller wavelengths, all 


smalle'h wavelengths should have reached the nonlinear limit too. 


For a rising bubble o r jet of plasma of form approximated by Eq.. (5) 


the wavelength A = R/& has an amplitude of approximately 1/2 s o that during 


deceleration by the field by Eq. (26) one generation further. growth can take 


place within the nonlinear limit. Therefore, the plasma should penetrate 


the field with a mixing length R/n. 


The subsequent motion of the plasma will be an expansion along the 


lines of force, since in the parallel direction less stress is exe,rted by the 


field. 
The principle direction of this expansion will be away from the hemi- 


sphere of detonation because .of tke following stresses, all of which tend to 


direct the plasma from a northern hemisphere explosion towards a southern 


hemisphere deposition: 
, 
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1. The verticle motion of initial expansion Eq. (5) has a component 


of. momentum along. the field lines towards the far conjugate point. 


2. 
The ambient field initially diverges in a direction towards the far 


conjugate point. 


3. The distorted field diverges more strongly.because of the ' line 


tying' in the partially ionized atmosphere at the near conjugate 


PO int . 


For these reasons the plasma will tr'avel along the lines of force until they 


intersect the earth's atmosphere at the far conjugate point with a distribution 


determined by the magnetic flux surfaces intersection with the earth. Since 


we have already shown that due to the low density of the plasma a t the time 


of mixing the interaction should be energy-cons erving , the velocity along the 


lines will be sound speed, namely, u 
. 
ex 


DEPOSITION IN THE ATMOSPHERE AT THE CONJUGATE POINT 


On intersection with the atmosphere at the conjn~a.t.e pnint, the plasma 


will have a width (longitude) a s weli a s latitude determined by both the mixing 


length R/T and the diamagnetic pressure. Since the plasma expa,nds intn a 


tube roughly 6000 k m long, the radius for a pressure the same as the bubble 


of R = 1000 km is also approximately 300 km, The altitude of deposition will 


be determined by the electron dynamic friction and the niicl.car Coulomb 


scattering. The impact velocity u 
determines which of these two dynamic 
ex 
frictions will dominate. If the atmosphere is not heated by impact to a 


temperature high enough s o that the free electron density is comparable to the 


nuclear density, then the electron dynamic friction becomes small for an 


impact velocity such that the relative kinetic energy of an electron is less 


than the binding energy in the atom, namely: 
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Below this .velocity the bound-free electron dynamic friction becomes 


small and so.the heating of the atmosphere is small until the nuclear scattering 


becomes .large. By then .the incident plasma is .stopped and the question of 


heating.the atmosphere to ionization becomes unimportant. 
Only in,the 
. 


8 
velocity range u 
= 1.5 x 10 does the question of heating become important. 
ex 


The energy available is reduced to W /3 due to radiation loss during. ioniza- 
t 
2 
tion and recombination, and this energy is deposited over an area ( ~ / 3 ) . 


Assuming twice the ionization potential in electron volts is required to ionize 


each atom because of radiation loss, the maximum density N max that could 


be ionized is : 


h N 
2 E . k = wt/3(R/3) 2 
max 
1 


for 


E. = 7 ev. 
1 


N 
9 
= 1.2 X 10 /cc. 
max 


This is too low a density to stop the plasma by electron free-free 


dynamic friction by Eq. ( l l ) , s o the dominant stopping,will be by bound-free 


8 
8 
for u 
> 1.5.X 10 cm/sec and by nuclear-coulomb for u 
< 1.5 .X 10 . 
ex 
ex 


The nuclear-Coulomb dynamic friction is given by Eq. (12) on the basis that 


both nuclei of a collision a r e unscreened by bound electrons. 
This is the 


case if the impact parameter ( ~ r / ~ ) ' / ~ 
is small compared to the orbit of the 


k-shell electron. 
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T h e minimum incident energy for which Eq. (12) is valid is determined 
. 
. 


by 


where 
Ry = Rydberg 
. 
. 
constant = 13 ev 


and 
i n n = 10 


giving 


L 
M u 
/2 = 13 kev. 
exp 


At lower incident energies the t e r m that sums the multiple scattering 


(anA) approaches unity and the cross section approaches the classical 
2 
Rutherford cross section a t M u 
/2 = 4 kev. 
Therefore an approximately 
ex 
-17 
2 
constant stopping cross section of u = 1.5X 10 
c m can be given for the 


L 
energy range 4 < - Mu /2 < - 13 kev. 
This corresponds to deposition at an 


altitude of 250 km. 


Since for the very much higher energy region the bound-free dynamic 
, 


friction should stop the incident plasma a t the same altitude (Eq. (11); 


u 
. = u /2), we can feel reasonably confident of a deposition height of 250 
exp 
0 


km. The distribution .according .to magnetic flux surface should therefore be 


compared to the earth's field at this altitude. 


MAGNETIC FIELD GR.ADIENT DRIFT 


In addition to .the spreading in flux surface due to the Taylor .instability 


mixing, there will also be a drift in longitude due to finite Larmor orbit drift 


in the gradient of the magnetic field. This depends upon the following 


parameters and assumptions: 


i 
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1. 
The region of most effective field gradient may .be the small' region 


.of high gradient in .the initial field distortion or the dipole gradient effective 


during the transit to the conjugate point. 


.2. The transverse electric field established by charge separation may, 


o r may not dominate the O B drift. 


3. The transverse ion velocity and equilibrium charge state effect both 


the' 
X 
and VB_ drift. 


Unfortunately, both regions of field - the distortion loop and dipole 


gradient - a r e effective in producing drift with the greatest uncertainty con- 


cerning the detailed loop shape and gradients. In addition the relative im- 


portance of the electron stress versus ion s t r e s s in the formation of the V B 2 


of the loop is dependent upon a detailed knowledge of the nonlineai.,instability 


mixing behavior. If the interaction boundary is a stable IM' layer as al- 


ready discussed, then .the ion stress is balanced by an electric field from 


charge separation and only the electrons drift transverse .to.the magnetic 


field. 
On the other hand, if the layer is unstable, the stress will be divided. 


according,to the respective temperatures -which in turn a r e dependent upon 


the form of instability. 


The drift in the dipole field and in the distortion loop will be calculated 


on the a priori assumptions: 


I 


' 
(1) The electron temperature is determined by the linear processes ': 


- 


only and will be taken as, 10 ev from Eq. (24) for the limited initial density 


(2) The ion trar~sver 
s e ,velocity i i ~ 
the earth's ficld is the expansion 


velocity u ex 
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F o r these assumptions the drift velocity v 
becomes 
D 


r 
VB/B 
V~ = Uex L 


r~ = ion Larmor radius = M u 
/ Z B. 
ex 
e 


The charge state Z is determined by the recombination rate in the 


expanded plasma a t the density n ( h / ~ ) ~ 
and Te 1 10 ev. 
This rate (for 
0 


3-electron configurations) is given by Post [post, 19611 a s 


-16 4 
0 8 X ; l O 
Z 


nuv = n 
e ff 


3/2 
k3 (Y 


e 


where k (y) will be close to' unity for kT' << E., the ionization potential of 
3 
e 
I . 


the state Zeff. 


For a recombination time of 1 sec, 


for T;e = 10 eu 


g i v h g 
z e f f % 1 i~ 
i~ 
. 


This means that for the density range considered, radiative recombination 


is much too slow a process to occur within the available transit time of the 


plasma. 
Consequently, the charge state will be determined by the maximum 


charge stripped state during the initial detonation products -air interaction. 


This in turn depends ,upon the runaway electron temperature in this region - 


and therefore upon the form of the interaction. From Eqs. (9) and (10) we 


will assume the a i r has been completely stripped to a mean charge of 7. The 


subseque~lt: charge exchange with neutral air atoms in traversing a distance 


4 
of the. earth's radius at a n ambient density of 10 /cc is small. F r o m Mott 


and Massey [ ~ o t t ' a n d 
Massey, 19491 


u = 10 
8 2 
8 
- I 6 (U 
/2.5 X 10 ) where u 
< 2.5 X 10 cm/sec. 
ex 
ex - 
(32) 
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The,fraction charge exchange to Z = 6 is less than 


-4 
,f < nu L = 6 X 10 
- 


8 
( L = earth's radius = 6 X 10 cm) 


and so is negligible. 


The Larmor radius then becomes 
- 4 
- 4 
r L = 5 . 3 X 1 0 / B U 
= 7 X 10 
u 
c m i n e a r t h ' s f i e l d . 
ex 
ex 


To calculate the drift in the earth's field the dipole field gradient is 


OB/B = 3 
/ 
~ 


so that the drift velocity becomes: 


and the drift D becomes 


-3 
D = L 
~ 
~ 
/ 
u 
~ 
~ 
= 2 x 10 
u 
cm. 
ex 
(33) 


Some fraction of the detonation products a s a precursor .may not have been 


subjected to the runaway temperature region ,of the mean interaction and so 


may have a much lower state of ionization. Choosing the extreme of charge 


one for a fission product the drift becomes 
- 1 
D = 10 
u 
cm. 
ex 


This distance is smaller than the mixing length spread provided u 
< u 
ex 
0 
and so should not be observable.. 


The drift in the formation of the distortion loop may be considerabdy 


larger. Again assuming the stress on the field is created by the ions, then 


the logarithmic gradient of the field should be roughly the ion Larmor radius 


for a length of time corresponding.to ,the intermixing of plasma and field. 


Since the drift velocity becomes roughly half the incident velocity uex9 and 


the period of drift is R/T u 
the total drift becomes 
ex' 
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BETA INJECTION FROM RADIOACTIVE DEBRIS 


The injection of energetic electrons from the p decay of the radio- 


active debris should occur during,the transit of the debris to the conjugate 


point. 
Only the decay after the unstable mixing of the plasma and field should 


permit injection .onto flux surfaces other than .the one intersecting,the detona- 


/' 
tion point. 
Since the time to reach the unstable mixing condition R/U 
is in 
ex 
general shorter (< 1 sec) than the major fraction of the p decay (- 3 .set), and 


the flight time to the conjugate point 1nn'gc.r than 5 secondni. thcn IIIC,.S~ 01 L21e 


decay p's will be injected into flux surfaces 'above the atmosphere. However, 


only those flux surfaces that do not intersect the earth's atmosphere for a 


reasonable mir'ror ratio will show long p lifetime. This minimum L surface 


[ ~ c n w a i n , 196 1 ] is approximat'ely L = 1.3 s o that only that fraction of the 


debris reaching this L surface and above:by unstable mixing should contribute 


to the long-lifetime P injection. 


NEUTRAL AND'SINGLE-CHARGED DETONATION PRODUCT BEHAVIOR 


. -. 
The considerations so far have.been for an average behavior and only 


the dominant interactions have been. considered. 
It is evident, however, that 


the charge state will fluctuate and the effects of a statistically small fraction 


may be important. In particular, the singly-charged high energy tail of the 


detonation products may drive the magnetic field and partially (X-ray) ionized 


air away from the remaining neutral fraction before ionization is complete, 


' and secondly, a small fraction of the debris neutralized by cha.rge exchange 


will expand across the field lines to give p-decay injection a t higher L values 


than would be reached by the Taylor instability mixing. 
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BETA INJECTION BY NEUTRAL DEBRIS 


.", 


The fraction of the debris that escapes a s neutral atoms should be 


generated predominantly by the slower .moving debris that expands upward 


after intersecting-with the air. The higher velocity precursor debris has a 


smaller probability for escaping as neutrals. The cross section for .pickup 


. from neutral a i r is [ ~ o t t 
and Massey, 19491 


= 10 -15 
2 
8 
cm f o r u 
< 1 0 . 
ex - 


On the other hand, the a i r through which the debris is rising and expanding 


will be ionized by the debris. Only that density region which is not fully 


ionized by the expanding debris will contribute to the charge neutr.alization. 


Assurning,the debris 'temperature during this phase of expansion is greater 


than.5 ev, the ionization rate will be temperature-independent and the ex- 


panded density determines when the further ionization of ambient neutrals 


is unlikely. 
This density is determined by 


-. 


giving 


The expansion of the plasma a t initial density 6 n to the density n 
0 
e 
takes place in a distance r into a solid angle of 1 radian such that: 


and the ambient neutral density becomes: 


T, 


N 
= n exp [-(l 
n /u 
n 
0 
0 
ex 1 lI3] . 
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he lieutral fraction becomes 


This is a major fraction only for n 
e lo8, o r an sltitudc greater than 550 km. 
0 


The p injection from this neutralized debris will occur according to 
- 1.2 
the usual t 
power law for fission product decay. It should be noted that 


for fairly'low altitude detonations where u 
is very small.(Eq. (20)) - of the 
. . 
C3L 


6 
order 10 cm/s ec'- neutralization will take place by recombination within 


the expanding bubble. 
This debris will be rising slowly s o that both p decay 


and sunlight will cause reionization. The decay spectrum will change within 


this time to give a different injection spectrum a s a function of L value. In 


addition, a significant change in the trapped P spectrum should be observed 


for flux surfaces reached in times of the order of a n hour, and the p velocity 


distribution function a t the equator will depend upon the debris reaching 


their corresponding geopotential. 


AURORAL EFFECTS 


? 


.The deposition of the detonation products and air plasma at the far 


co.njugate point will heat . up 
. the a i r in which it i s deposited, causing the air 


to expand upward and intersecting some of the trapped p-decay electrons. 


The optical pattern will correspond to the illumination of the fluxsurfaces 


on which high-energy electrons a r e trapped and hence to ,the gradient field 


drift separated charge state of the @-decaying debris. 


As previously discussed, the heating takes place a t 250 km, depositing 


W 
electron- volts per atom during.the transit time of the plasma. 
This 
D 


energy is: 
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(4 1) 


10 
for ~ 
/ 
3 
= 300 km, n = 2.X 10 
and 


then WD = 40 ev/atom. 


The transit time is s o long, however, (up to.tens of seconds) that this energy 


will be almost entirely radiated away a s it is deposited. 
The residual 


temperature will therefore be determined by the radiation properties of the 


air - which becomes .transparent at roughly 1 /2-ev temperaturL. The 


vertical expansion velocity corresponding ,to twice sound speed at this. tempera- 


5 
ture will be approximately 3 X 10 cm/sec, permitting,the air to reach 500 ' 


to 1000 km. altitude in several minutes. The optical radiation corresponds 


to the excitation of the a i r with approximately 10% fluorescent efficiency for 


an energy deposition rate corresponding to the usual minimum ionization 


'2 
:rate of 2 .Mev /g/cm . 


The p-decay electron density for uniform phase-space injection into a 


magnetic tube of force 300 km wide and 6000 km long for an injection time 


longer than .5 seconds becomes: 


= W /E 
Vol = 2.5 X lo3 electron/cc 
P 
P 


where 


W. 
= p-decay energy = 0.05 W 
p 
. 
.:. 
t 9 


E = mean p energy = 1 Mev, and 
P 


Vol = Volume of flux tube. 


The radiated energy loss rate then becomes 


6 
-6 
ER = cn P X 2 X 10 k. = LO 
nrgs/cr. 
sec 
P 
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where c = electr.on velocity 


p = idensity after. 10-fold expansion 


= fluorescent efficiency ,- 10%. 


The total optical brightness Br f o r a transparent thickness of 300 km 


becomes 
7 
2 
B 
= 3.X 10 E 
= 30 erg/cm 
sec. 
r 
r 


The heating rate of the air is ten-fold bigger than Eq. (43) and corresponds 
. 
I 


to 1 ev per atom during .the 100 seconds of expansion. This is enough to , 


enhahce th{ expansion so that sdme air would, be expected to ionize and expand 


back along the lines of force from both conjugate regions.to meet at the 


equator. . The energy loss rate per electron is 2000 ev/sec, but during.the 


available time before the gradient field drift carries the electrons eastward 


the fractional energv loss, i s small; The energy lnss ncciirs f n r a ti-?: 


7 5 ( . R / ~ ) , / v ~ 
= 30 sec where 
. . 


6 
VD = 1 0 cm/scc from Eq. (27) for a 1-Mev electron 


go-that only the low energy electrons with larger energy loss rate and, slower 


drift will be lost in the rising air mass: 
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CP 
Fig. 1. The time sequence of the expansion of detonation products in a stratified atmosphere is 
represented by the approximate snowplow shock contour. The air and detonation products in this case 
P 
(D 


a r e considered neutral and only subs equent fission decay gives debris interacting with magnetic field. 
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Fig. 2. 
Ionization magnetic shock (moving frame). 


MAGNETIC 


Fl ELD INTO PAGE 


LEGAL NOTICE 


This report was prepared a s a n account of Government sponsored work. 
Neither the United States, nor the Commission, nor a n y person acting on 
behalf of the Commissiori: 


A. Makes any warranty or representation, expressed or implied, with 
respect to the accuracy, completeness, or usefulness of the information con- 
tained in this report, or that the use of any information, apparatus, method, 
or process disclosed in this report may not infringe privately owned rights; or 


B. Assumes any liabilities with respect t o the u s e of, or for damages 
resulting from the use of any information, apparatus, method or process dis- 
closed in this report. 


A s used in the above, "person acting on behalf of the Commission 
includes any employee or contractor of the commission, or employee of such 
contractor, t o the extent that s u c h employee or contractor of the Commission, 
or employee of s u c h contractor prepares, disseminates, or provides a c c e s s 
to, any information pursuant t o his employment or contract with the Commis- 
sion, or his employment with s u c h contractor. 


